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Abstract 

^^N(p, 7)^^0 is one of the key reactions in the hot CNO cycle which occurs at stellar tempera- 
tures around Tg > 0.1. Up to now, some uncertainties still exist for the direct capture component 
in this reaction, thus an independent measurement is of importance. In present work, the angular 
distribution of the ^^N((i, n)^^0 reaction at Ec.m. = 8.9 MeV has been measured in inverse kine- 
matics, for the first time. Based on the distorted wave Born approximation (DWBA) analysis, the 
nuclear asymptotic normalization coefficient (ANC), C'i\^2' ground state of ^^O ""^^N + 

p is derived to be 5.42 it 0.48 imT^^"^. The ^^N(p, 7)"'^^0 reaction is analyzed with the R-matrix 
approach, its astrophysical S-factors and reaction rates at energies of astrophysical relevance are 
then determined with the ANC. The implications of the present reaction rates on the evolution of 
novae are then discussed with the reaction network calculations. 
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I. INTRODUCTION 



In stellar-evolution models, hydrogen burning in massive stars proceeds largely through 
the CNO cycle. For the normal CNO cycle, the dominant sequence of reactions is 
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When temperature increases, the decay of ^^N limits the cycle, and most of the C, N 
and O nuclei would be processed into ^^N. Consequently, the ^^N(p, 7)^^0 reaction provides 
a second channel for destruction of ^'^N, and the dominant sequence becomes 



This reaction sequence is called hot or "/5-limited" CNO cycle, and the decays of ^''O and 
^^O limit this cycle. The CNO cycles convert four hydrogen nuclei into an alpha particle 
and the energy release in the cycles is about 26.7 MeV, which is the important source of 
stellar energy generation [l]. Since the decays of ^^O and ^^O are much quicker than 
that of ^^N, the hot CNO cycle should produce energy much faster than the normal CNO 
cycle. Hence, a rapid change of the temperature dependent energy generation rate occurs 
when the CNO cycle transits from the normal one to the hot one. ^^N(p, 7)^^0 is one of the 
important reactions which controls this transition [S]. Therefore precise determination of the 
rates for the ^^N proton capture reaction is vital for determining the transition temperature 
and density between the normal and hot CNO cycles. 

At the energies of astrophysical interest, the ^^N(p, 7)^^0 reaction is dominated by the 
low energy tail of the s-wave capture on the broad 1~ resonance at Er = 527.9 keV (which 
has a total width of 37.3 ± 0.9 keV). A considerable effort has been expended in recent 
years to determine the parameters for the resonance. These include the direct rneasure 
ments using the radioactive ^^N beam j^, 4], particle transfer reactions [s, ^ 7, 
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, and 
10, 



111. 



Coulomb dissociation of high-energy ^^O beams in the field of a heavy nucleus 
The direct capture contribution is significantly smaller than the contibution due to the tail 
of the resonance within the Gamow window. But since both resonant and non-resonant 
captures proceed via s-waves and then decay by El transitions, there is an interference 
between the two components. Thus the capture reaction within the Gamow window can 
be enhanced through constructive interference or reduced through destructive interference. 
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The non-resonant component of the cross section has been calculated by several groups, 
either separately or as part of the calculation of the total cross section l|, Q, Q, Q- Since 
there are significant differences among the various calculations, the determination of the 
^^N(p, 7)^^0 direct capture component through an independent approach is greatly needed. 
A practicable method is to extract the direct capture cross section of the ^^N(p, 7)^''0 reac- 
tion using the direct capture model Q and the spectroscopic factor (or ANC), which 
can be deduced from the angular distribution of one proton transfer reaction. Decrock et 
al. extracted the spectroscopic factor for ^^O ^^N + p from the ^^N((i, n)^^0 cross section 
0. Tang et al. derived the ANC for ^^O ^ ^^N + p from the ^^N(13n,^^0)13C angular 
distribution [18]. The S-factors for the direct capture of the ^^N(p, 7)^^0 reaction from 
these two works differ from each other by a factor of 30%. Thus, further measurement is 
important for the determination of the spectroscopic factor (or ANC) for ^^O — > ^^N + p 
and the astrophysical S-factor of the ^^N(p, 7)^^0 reaction. 

In the present work, we have measured the angular distribution of the ^^N((i, n)^^0 re- 
action at Ec.m. = 8.9 MeV in inverse kinematics. The spectroscopic factor and ANC were 
derived based on distorted wave Born approximation (DWBA) analysis, and used to calcu- 
late the astrophysical S-factors and rates of ^^N(p, 7)^^0 direct capture reaction at energies 
of astrophysical interest with the R-matrix approach. We have also computed the contri- 
bution from the resonant capture and the interference effect between resonant and direct 
capture. The total reaction rates are then used in the reaction network calculations at the 
typical density and temperature of novae environment. 



II. MEASUREMENT OF THE ^^^((1,71^^0 ANGULAR DISTRIBUTION 

The experiment was carried out using the secondary beam facility Q, Q of the HI- 13 
tandem accelerator, Beijing. An 84 MeV ^^C primary beam from the tandem impinged on a 
4.8 cm long deuterium gas cell at a pressure of 1.6 atm. The front and rear windows of the 
gas cell are Havar foils, each in thickness of 1.9 mg/cm^. The ^^N ions were produced via 
the ^H(^^C, ^^N)n reaction. After the magnetic separation and focus with a dipole and a 
quadruple doublet, the secondary beam was further purified with a wien filter. The 69 MeV 
^^N secondary beam was then delivered with typical purity of 92%. The main contaminants 
were ^^C ions out of Rutherford scattering of the primary beam in the gas cell windows as 
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FIG. 1: Schematic layout of the experimental setup 



well as on the beam tube. The ^^N beam was collimated with two apertures in diameter of 3 
mm and directed onto a (CD2)„ target in thickness of 1.5 mg/cm^ to study the ^H(^^N/^0)n 
reaction. The typical beam intensity and beam energy spread on the target were 1500 pps 
and 1.8 MeV FWHM for long-term measurement, respectively. A carbon target in thickness 
of 1.5 mg/cm^ served as the background measurement. 

The experimental setup is shown in Fig. [H A 300 /im thick Multi-Ring Semiconductor 
Detector (MRSD) with center hole was used as a residue energy {Ej.) detector which com- 
posed a AE — Ej- counter telescope together with a 21.6 fim thick silicon AE detector and a 
300 yum thick silicon center Er detector. Such a detector configuration covered the laboratory 
angular range from 0° to 5.4°, and the corresponding angular range in the center-of-mass 
frame is from 0° to 66.5°. This setup also facilitates to precisely determine the accumulated 
quantity of incident ^^N because the ^^N themselves are recorded by the counter telescope 
simultaneously. 

The accumulated quantity of incident ^^N is approximately 3.54 x 10^ for the (CZ^2)n 
target measurement, and 1.18 x 10^ for background measurement with the carbon target. 
Fig. [2] (a) - (d) display the AE — E^ scatter plots for the first four rings, respectively. For 
the sake of saving CPU time in dealing with the experimental data, we set a cut line of AE 
= 19 MeV. All the events below the line are scaled down by a factor of 1000, and the ^^O 
events are not affected by this cut. The four two-dimension gates plotted in Fig. [2] (a) - 
(d) are the ^^O kinematics regions based on the Monte Carlo simulation, taking the beam 
spot size, energy spread, angular divergence and the target thickness into account. The ^"^O 
events can be clearly identified through this figure. Fig. [3] displays the comparison of the 
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FIG. 2: (Color online) Scatter plots oi AE vs. for the ^^N{d,nY^O reaction measurement with 
{CD2)n target, (a) - (d) display the IS.E- E^ spectra for the first four rings of MRSD. 

events from (CD2)n target with the background from carbon target in the ^^O kinematics 
regions for the first four rings. The background events in the first ring of MRSD mainly 
come from the pileup of ^^C contaminants in the beam, they disappear in the outer rings. 
After the background subtraction, the angular distribution in center of mass frame for the 
forward angles is given in Fig. |H The uncertainties of differential cross section mainly arise 
from the statistics, the assignment of ^^O kinematics regions, the uncertainties of the target 
thickness and the solid angle. The angular uncertainties include the random reaction point 
in the target, the angular uncertainty of the secondary beam, the angular straggling of ^^N 
and ^^O in the target and the AE detector. The total angular error for each ring is about 
0.6 degree, less than the width of each ring. 



III. DETERMINATION OF THE NUCLEAR ANC 

The spins and parities of ^^N and ^^O (ground state) are 1/2" and 0^, respectively. The 
cross section of the ^^N((i, n)^^0 reaction is dominated by the s-wave proton transition to 
lpl/2 orbit in ^^O ground state. If the reaction is peripheral, the differential cross section 
can be expressed as 

da /C'rfN2/ 1,1/2x2 (n\ f-W 

a\l Od ^1,1/2 
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FIG. 3: (Color online) Comparison of total energy spectra between (CD2)n and pure Carbon target. 
(a)-(d) represent the spectra for ring 1,2,3 and 4. The solid bar and the empty bar stand for the 
^^O spectra from (CD2)n target and carbon target, respectively. 
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FIG. 4: The angular distribution of the '^^N{d,nY'^0 reaction at Ec 
DWBA calculations using different optical potential parameters. 



8.9 MeV, together with 



where {^)exp and crij{9) denote the measured and theoretical differential cross sections 



respectively. stand for the nuclear ANCs for the ^^O ^^N + p and d ^ n 

+ p virtual decays, fe"{/2 ^^'-^ being the single particle ANCs of the bound state protons 
in ^^O and deuteron. By knowing the value of Cd , the Cl*y2 '^^^ then be extracted by 
normalizing the theoretical differential cross sections to the experimental data by Eq. ([T]). 

The DWBA code DWUCK [21] is used to compute the angular distribution. All the op- 
tical potential parameters for the entrance channel are taken from Ref . 2^ , the ones for the 
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FIG. 5: Comparison of spectroscopic factors with the ANCs derived from the present experiment 
for different geometries of the Woods-Saxon potentials. 



exit channel are from Refs. j22j] and 23|], respectively, these parameters are listed in Tabled 



In the present DWBA calculation, the differential cross sections at three forward angles are 
used to extract the ANC, and Cd is taken to be 0.872 fm~^/^ from Ref. 2^. The normalized 
angular distributions from the six sets of optical potential parameters are also presented in 
Fig. m each curve corresponds to one nuclear ANC, C^i^j^, the spectroscopic factor is calcu- 
lated with The nuclear ANCs and the spectroscopic factors deduced from the present 
experimental data are listed in Table HTl the average values of them are 5.42 ± 0.48 fm~^/^ 
and 1.88 ± 0.34, respectively. The present ANC accords with the result extracted from the 
^^N(^^N, ^'^OY^C transfer reaction by Tang et al. [18], and the present spectroscopic factor 
is larger than the previous one (0.9) extracted from the total cross section of ^'^N((i, n)^^0 
at lower energy V]\. The uncertainties of the nuclear ANC and the spectroscopic factor 
are mainly from the difference of the calculated angular distributions with different optical 
potentials, as well as the experimental errors. Since we do not measure the optical potential 
parameters and used six sets of potential parameters from the neighboring nuclei, the error 
bar of present work is a bit larger than that of Ref. 18| . Fig. [5] shows the comparison of the 
spectroscopic factors with ANCs of ^^O ^'^N + p from the different geometry parameters 
of the Woods-Saxon potential for the single particle bound state (by changing the radius 
and diffuseness tq and a). One can see that the spectroscopic factors vary strikingly, while 
the ANCs are nearly a constant, thus indicating that the ^^N((i, n)^^0 reaction at present 
energy is dominated by peripheral process. 
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TABLE I: Optical potential parameters used in DWBA calculations, where V , W are in MeV, r 
and a are in fm, the geometrical parameters of single particle bound state are set to be tq = 1.25 
fm and a = 0.65 fm. Dl, D2 and D3 correspond to the optical potentials for d + ^^N, and Nl, N2 
represent the ones for n + ^^O. 



Channel 


Dl 


Entrance 
D2 


D3 


Nl 


Exit 

N2 


Vr 


117.9 


116.0 


130.4 


49.2 


61.56 


ror 


0.81 


1.0 


0.9 


1.2 


1.14 




1.07 


0.8 


0.9 


0.65 


0.57 


Wv 




4.13 












1.0 








Clujv 




0.8 








Ws 


19.61 


4.13 


6.63 


6.0 


7.74 


ros 


1.84 


2.0 


1.90 


1.2 


1.14 


as 


0.35 


0.6 


0.56 


0.47 


0.5 


Vso 




6.76 




7.0 


5.5 


roso 




1.0 




1.20 


1.14 


aso 




0.8 




0.65 


0.8 


roc 


0.81 


1.5 


1.30 







IV. ASTROPHYSICAL S-FACTOR OF i3N(p, 7)1^0 



According to the traditional direct capture model jl5, 16|, l25|, the direct capture of the 
^^N(p, 7)^^0 reaction is believed to be dominated by the El transition from incoming s wave 
to bound p state. The direct capture cross section can be computed by 

100 = —k^^e'^A^SiJ J rfrrV,(r)^Jr)|2. (2) 

where = e^/hc is the wave number of the emitted 7-ray (of energy e^) and e = eN/A 
is the El effective charge for protons, Aij corresponds to the angular part depending on 
the initial and final angular momenta of the transition, Sij is the spectroscopic factor of 
the configuration ^^O + p, ^pi^{r) is the bound state wave function of the relative 
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TABLE II: The O nuclear ANC, ^j^^ and spectroscopic factor, ^^g; deduced from the angular 
distribution of the ^^N(d, n)^^0 reaction using the combination of optical potentials for the entrance 
and exit channels. 



optical 
Dotpntials 


1,1/2 


c-"0 
1,1/2 


Dl-Nl 


h.Tl ± 0.42 


1.77 ± 0.28 


D1-N2 


4.95 ± 0.17 


1.56 ± 0.11 


D2-N1 


6.02 ± 0.61 


2.31 ± 0.47 


D2-N2 


5.42 ± 0.29 


1.87 ± 0.20 


D3-N1 


5.56 ± 0.27 


1.97 ± 0.19 


D3-N2 


5.31 ± 0.19 


1.80 ± 0.13 


average 


5.42 ± 0.48 


1.88 ± 0.34 



5 I 1 1 1 1 1 r 




■2 I 1 I 1 I 1 I 1 

5 10 15 20 

Radius (fm) 



FIG. 6: The integrand of the E\ transition matrix element based on a single-particle model at 
resonant energy. 

motion of p+^^N in ^^O calculated in the Woods-Saxon potential, ipi^f^) is the optical model 
scattering wave function of the colliding proton and ^^N. If the spectroscopic factor 5*;^- is 
deduced from the ^^N{d,nY^O transfer reaction, the ^^N(p, 7)^^0 cross section can then be 
calculated by Eq. ([2]). 

However, this is not the case here, as a result of the tight binding of the last proton in ^^O, 
the contribution to the ^^N(p, 7)^^0 direct capture reaction at small r in Eq. ([2]) is important. 
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The integrand of the El transition matrix element at resonant energy is calculated based 
on a single-particle model, as shown in Fig. [61 One can see that the contribution at small r 
is of significance, the simple direct capture model may be not valid due to the many particle 
effects. In this case, the integral is very sensitive to the optical potential parameters and the 
spectroscopic factor required for Eq. ([2]) has significant uncertainties, as can be seen from 
Fig.O 

In this work, we will use R-matrix method to avoid the above problems. For the radiative 
capture reaction B + b A + 'y, the R-matrix radiative capture cross section to a state of 
nucleus A with a given spin Jf may be written as [26 1 



^ TT 2Jj + 1 2 /A\ 

Ili 

Here Jj is the total angular momentum of the colliding nuclei B and b in the initial state, 
and Jb are the spins of nuclei b and B, and /, k, and /j are their channel spin, wave number 
and orbital angular momentum in the initial state. Uu.jjji is the transition amplitude 
from the initial continuum state {Ji,I,li) to the final bound state {Jf,I). In the one-level, 
one-channel approximation, the resonant amplitude for the capture into the resonance with 
energy Er^ and spin Jj, and subsequent decay into the bound state with the spin Jf can be 
expressed as 

' E-Er^+i^ 
Here we assume that the boundary parameter is equal to the shift function at resonance 
energy and 0;. is the hard-sphere phase shift in the /jth partial wave. 



arctan 



(6) 



where Ff, and Gf^ are the regular and irregular Coulomb functions, is the channel radius. 
The Coulomb phase factor uoi- is given by 

h 

ooi, = ^arctan(— ), (7) 

n=l ^ 

where r/i is the Sommerfeld parameter, rjfj;. (i?) is the observable partial width of the res- 
onance in the channel B + b, r'^'j^{E) is the observable radiative width for the decay of 
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the given resonance into the bound state with the spin Jf, and Tj. Yli^biu 
servable total width of the resonance level. The energy dependence of the partial widths is 
determined by 



^ blli 



PiAE) 



^tnXERn) 



and 



_ E + Ef . 2L+lpJi I 



XE) 



ErJ, 



(8) 



(9) 



where F^}^. and r'^j^{Eii^) are the experimental partial and radiative widths, £/ is the 
proton binding energy of the bound state in nucleus A, and L is the multipolarity of the 
gamma transition. The penetrabihty Pl^ (E) is expressed as 

kVr 



The nonresonant amplitude can be calculated by 



(10) 



TjNR 



_^2)3/2,;i»+L-i/+lpi(i-'ii 



1 L+1/2 



hk 



X 



mt 



mo 



(L + 1)(2L + 1) 



xGi^ik, re) Wi^ {2Kr,) ^^(/i OLO \lfO) 
xUiLlfJJ;hJf)J'L{klf), 



(11) 



where 



— f 

c Jr. 



GiXKr,) 



WipnT,)VFiXk,T,) 



(12) 



Here, Wi{2nr) is the Whittaker hypergeometric function, k — ^2\iBh^f and // are the wave 
number and relative orbital angular momentum of the bound state, and k^ = {E + ef)/hc 
is the wave number of the emitted photon. 

The non-resonant amplitude contains the radial integral ranging only from the channel 
radius Tc to infinity since the internal contribution is contained within the resonant part. 
Furthermore, the R-matrix boundary condition at the channel radius Tc implies that the 
scattering of particles in the initial state is given by the hard sphere phase. Hence, the 
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FIG. 7: Astrophysical S-factors as a function of Ec.m. for the ^^N(p, 7)^^0 reaction. The dotted 
hne is the contributions from the direct proton capture. The sohd, dashed, and dashed-dotted Unes 
indicate the total S-factors from the present work, Refs. [l^ and [l7], respectively. 

problems related to the interior contribution and the choice of incident channel optical 
parameters do not occur. Therefore, the direct capture cross section only depends on the 
ANC and the channel radius r^. 

The astrophysical S-factor is related to the cross section by 



S{E) = Ea{E) expiEc/E) 



1/2 



(13) 



where the Gamow energy Eq = 0.978Z^Z|/i MeV, /i is the reduced mass of the system. 
According to the experimental ANC (5.42 ± 0.48 fm~^/^) from the present work, and the 
resonance parameters {Eji = 527.9 ± 1.7 keV, TtotiEu) = 37.3 ± 0.9 keV, and T^^En) = 
3.36 ± 0.72 eV) from Ref. the S-factors for direct and resonant captures can be then 
derived, as demonstrated in Fig. [71 

Since the incoming angular momentum (s-wave) and the multipolarity (£"1) of the direct 
and resonant capture 7-radiation are identical, there is an interference between the direct 
and the resonant captures. In this case, the total S-factor is calculated with (l5|] 



Stot{E) = Sdc{E) + Sres{E) ± 2[Sdc{E)Sres{E)Y^' COs(5), 

where S is the resonance phase shift, which can be given by 



(14) 



arctan 



2{E-Er)- 



(15) 
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TABLE III: The present total reaction rate for 


13n(p, 7)1^0, Na < 


(TV > (cm^mole ^s as a 


function of temperature, together with the previous results. 






Present work 


Kef. [langj 


NAHRF, 


0.01 


AAA X 10 


4.18 X 10 




0.02 


6.02 X 10 ° 


5.72 X 10 ° 


9 78 X 10~16 

^ • i O A XL/ 


0.03 


5.60 X 10 ^ 


5.35 X 10 ^ 


9 63 X 10~13 


0.04 


4.16 X 10 


3.98 X 10 


1 QQ X 10""ll 

± .ijij A XL/ 


0.05 


8.89 X 10 ^" 


8.53 X 10 ^" 


4:.Ot: a XL/ 


0.06 


9.19 X 10 ^ 


8.84 X 10 ^ 


t:.00 a XL/ 


0.07 


5.94 X 10 ° 


5.72 X 10 ° 


Q 09 V 1 

O.L/^ A XL/ 


0.08 


2.77 X 10 ' 


2.67 X 10 ' 


X.^^ A XL/ 


0.09 


1.02 X 10 ° 


9.86 X 10 ' 


^ 4^ V iri"''' 

fJ.^O A XL/ 


0.1 


3.15 X 10 ° 


3.04 X 10 ° 


1 71 V in~6 

X. 1 X A XL/ 


0.13 


4.41 X 10~^ 


4.27 X 10"^ 


Zi.fJU A XL/ 


0.17 


5.32 X 10 


5.16 X 10~^ 


^ ^4 V in~^ 

A XL/ 


0.21 


3.34 X 10 ^ 


3.24 X 10 


Z.ZZ A XL/ 


0.25 


1.44 X 10 


1.39 X 10 


Q V in~^ 

c/.0(J A XL/ 


0.29 


5.00 X 10~^ 


4.84 X 10~^ 


*^ 47 V 1 

0.'4:/ A XL/ 


0.33 


1.56 X 10~^ 


1.51 X 10~^ 


1 1 1 V in~i 

X . X X A XL/ 


0.37 


4.56 X 10 


A A -t f\ — 1 

4.41 X 10 


41 V in~i 

O.^X A XL/ 


0.41 


1.24 X 10" 


1.20 X 10" 


Q Q1 V in~i 

t/.C/X A XL/ 


0.45 


3.07 X 10" 


2.98 X 10" 


9 fifi V 1 n'^ 

Zi.UU A XL/ 


0.49 


6.87 X 10" 


6.69 X 10" 


fi fiQ V 1 n'^ 


0.53 


1.39 X 10^ 


1.36 X 10^ 


1 9S V If)! 

X.ZiO A XL/ 


0.57 


2.59 X 10 


2.54 X 10 


9 44 V 1 fl^ 

A XL/ 


0.61 


4.46 X 10 


4.38 X 10 


4 97 V 

^.Zi / A XL/ 


0.65 


7.20 X 10 


7.09 X 10^ 


fi QQ V 


0.69 


1.10 X 10 


1.09 X 10 


1 ns V in^ 

X.L/O A XL/ 


U. 1 o 


1 fin V ifi^ 


±.tJO A XL/ 


1 V Ifl^ 

X.fJO A XL/ 


0.77 


2.23 X 10^ 


2.22 X 10^ 


9 99 V in^ 

£j.£j£j A XL/ 


0.81 


3.01 X 10^ 


2.99 X 10^ 


3.01 X 10^ 


0.85 


3.94 X 10^ 


3.92 X 10^ 


3.95 X 10^ 


0.89 


5.02 X 10^ 


5.00 X 102 


5.04 X 10^ 


0.93 


6.26 X 10^ 


6.23 X 10^ 


6.30 X 10^ 


0.97 


7.64 X 10^ 


7.59 X 10^ 


7.70 X 10^ 
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Generally, the sign of the interference in Eq. (HM has to be determined experimentally. 
However, it is possible sometimes to infer this sign. The interference between the res- 
onant and direct capture contributions is constructive below the resonance energy and 
destructive above it, which has been observed from the isospin analog ^^C(p, 7)^^N* (2.31 
MeV) reaction \v)^. Recently, Tang et al. deduced constructive interference below the 
resonance using an R-matrix method [18]. Based on this interference pattern, the present 
total S-factor is then obtained. Fig. [7] shows the comparison of total S-factors from the 
present work, Refs. [l8| and [l7|. Our updated total S-factors are about 40% higher than 
the previous ones in Ref. l3] at low energies and is in good agreement with that in Ref. 18|. 



V. THE ASTROPHYSICAL REACTION RATE 

The astrophysical reaction rate of ^^N(p, 7)^^0 is calculated with 



xexp[-(|j)'/^-^]dB, (16) 



where is Avogadro constant. The updated rates are listed in Table UTTl together with the 
previous ones from Ref. [isl and NACRE's compilation. The results from the three works 
agree with each other within a factor of 2 at low temperature of T < 0.2 GK and are almost 
identical at higher temperature of T > 0.7 GK. 

The present total reaction rates as a function of temperature Tg (in unit of 10^ K) are 
fitted with an expression used in the astrophysical reaction rate library REACLIB 

Na{(Jv) = exp[-5.2635 + 0.0864X9-^ - 21.5656T9"^/^ + Se.OSTSTg^^ - 4.9432r9 + 0.3937T^^^ 

-9.7467 ln(r9)] + exp[108.6965 + 0.6657r9"^ - 47.9051T9~^/^ - 59.492ir9/^ + 5.OI45T9 
-0.2488X9/^ + 4.4288 ln(r9)]. (17) 

The fitting errors are less than 5% in the range from Tg = 0.01 to Tg = 10. 



14 



0.4 

FIG. 8: (Color online) Energy production rates of the CNO and hot CNO cycles at p = 500 (right) 
and 5000 (left) g/cva? for novae with the update ^^N(p, 7)^^0 reaction rates from present work and 
nacre's compilation. 

For a given density p, the reaction network equations and the energy source equation 
have the following forms: 

r,-F(r„T) = 

e + Y.NAM,^Yi = Q ^^^^ 

i 

where Yi are the nuclear abundances, e is the energy production rate per unit mass, i^j 
= 1,2,- ■ ■ ,N, and N is the number of nuclear species. F denotes nonlinear functions of 
the arguments, and MjC^ is the rest mass energy of species i in MeV. At equilibrium, the 
abundances do not change with the time approximately, i.e., Yi ~ 0, the energy production 
rate can then be calculated by substituting the reaction rates into Eq. (ITSl) . Fig. [8] shows 
the energy productions of CNO and hot CNO cycles at density p = 500 and 5000 g/cm'^ for 
novae with the ^^N(p, 7)^^0 reaction rates from present work and NACRE's compilation. 
One can see that the hot CNO cycle would begin to run earlier and produce more energy 
with our updated ^^N(p, 7)^^0 reaction rates. The present result shows that about 5% of 
additional energy could be produced at the temperature range from 0.07 to 0.15 GK, which 
implies that the evaluation of a novea may be affected. 
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VI. SUMMARY 



In summary, ^^N(p, 7)^'*0 is one of the key reactions which trigger the onset of the hot 
CNO cycle. We have measured the angular distribution of the ^'^N((i, nY^O reaction at -Ec.m. 
= 8.9 MeV, and deduced the nuclear ANC and spectroscopic factor for the ^''O ground state. 
The astrophysical S-factors and reaction rates for ^^N(p, 7)^^0 are then extracted with the 
R-matrix approach. Our result is in good agreement with that from the ^"^N (^^N/^0)^^C 
transfer reaction by Tang et al. |l8]. The reaction network calculations have been performed 
with the updated ^'^N(p, 7)^^0 reaction rates, the result shows that 5% additional energy 
could be generated through the CNO and hot CNO cycles at the typical densities and 
temperature range from 0.07 to 0.15 GK for the novae, this may affect the evaluation of 
novae. 
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